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Abstract 

We show that MS perturbation theory develops tachyonic singularities for some value 
of the dimensional regularization scale /x unless the physical Higgs mass exceeds some 
(cutoff dependent) minimum value. 

PACS: 14.80.Bn, ll.lO.Hi, 12.15.Lk 

Key words: Higgs, Boson, Mass, Bound, Lower 



* Physics Dept, University of Pittsburgh, Pittsburgh PA 15260; E-Mail:willey@vms. cis.pitt.edu 



1 



For the history of the subject we refer to reviews and quote only recent papers 
containing the latest refinements of the conventional approach [0,0. Once the experi- 
mental lower bound on the top-quark mass exceeded 80 GeV, attention shifted from the 
original Linde- Weinberg bound, based on the properties of the one-loop effective potential 
for small (p, close to the minimum, to the large (p behavior of the effective potential , as 
determined by renormalization group (RG) considerations. 
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Here, ^(t) is the anomalous dimension factor, and X{t) is the MS running coupling con- 
stant . 

It is then argued that vacuum stability requires A (t) > up to some high scale, ~ 
Mgut,otMpi, {Mq ~ Mz, or rrit, or 24QGeV).To implement this condition, one has to 
know (or approximate) the (3— function, integrate the RG differential equations starting 
from some initial values, A (0), (7^(0), and relate the smallest acceptable A (0) to a physical 
Higgs mass. (In the latest effective potential calculation Q the condition is Ae//(t) > 
where Ae// differs perturbatively from A ). 

We propose a new approach. It is also of the perturbative RG variety, but is not based 
on the effective potential . The results are found to be insensitive to scale ambiguity. The 
essential input is that one is perturbing about the correct vacuum. A neccessary condition 
for this is that the vev of the (shifted) field be zero, order by order in perturbation theory, 
and the MS renormalized mass squared in the MS propagator of the shifted field be 
positive. 

We start by computing the relation between the perturbative pole mass and the MS 
mass, for both the Higgs boson and the t-quark. The relation follows from the perturbative 
definition of the pole mass, 

= ReD'\M*') = M*'-M^-Re^{M*') (1) 

In this equation, D{q^)^ and S(g^) are the two- point Green Function and self-energy 
function, renormalized according to the MS prescription. M* is the perturbative pole 



mass. The result is 



M* =m'{1+A(3Joo(M* )+9J^(M* ))+iV,y2(_^_4— ^)/^(M* )+2(C,-l)} (2) 

m is the t-quark mass, and y is the t-quark Yukawa couphng (m = ■^). The contributions 
from the electroweak gauge sector, proportional to g2, di, have also been calculated, but 
are not written out here, they will be included below. The term Cv ~ ^ comes from a 
finite shift of the vev required in the MS scheme to enforce < H >= through one- loop 
order. It will cancel out of the ratio computed below, so we do not have to give its value 
here. is the dimensionally regularized MS scalar one-loop two-point integral. 

= 16^ [In ^ + /o^ dx In 

Then (0) is 



M*^ =m'{1 + -— [121n— -24+3v^7r]+Ar,-f^(l--)[ln^+/(r)] + 2(C-l)} (4) 



where 



M , 4-r2 / r2 

r = — , /(r) = — 2 + 2W — arctan - 



m ' ^ ' ' V V 4 — 

The corresponding calculation for the t-quark gives 
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l^[^ln^ + A(r)] + ^CH8-61n^' 



m*^ = m2{l + -|-^[- In — + A(r)] + ^C,.(8 - 6 In + 2(C„ - 1)} (5) 

where 



A(r) = -4 + — + (-r^ - -r"^) Inr^ -(4 - r^) 2 arctan \/ ^ J" 
y ' 2 ^2 4 ^ 2^ ' V r2 

We take the ratio of (0) to and expand to one-loop order. 

si = S{1 + Ti-[12 In f - 24 + 3y37r] + i[iV,(l - f + /(r)) - | In ^ - A(r)] 

+iCf (6 In ^ - 8) + ^72^, ^7? terms + 2 - /oop} 

(6) 

The Ct) ~ 1 terms, which also contain explicit dependence on In/x^^ have cancelled out. 



A necessary condition for the MS perturbation calculations to be defined in the broken 
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symmetry phase is that M'^^fn^ be positive. Since the ratio of pole masses is positive, 
satisfies the requirement perturbatively, for ^ around the weak scale. For large /i^, 
one has to provide a RG treatment of the large logarithms, just as in the conventional 
calculation involving the effective potential . In the broken symmetry phase, one can 
define the renormalized coupling constants such that the relation 

^ - 4 (T) 

is exact when all the quantities are either "star" (on-shell renormalization scheme) or 
"bar" (MS" renormalization scheme) [§| , P] . Thus, not all quantities in can be varied 
independently as functions of /i. we use (|^ to eliminate A appearing in (|]). To leading 
(one-loop) order, the scale dependence of the ratio of MS masses is determined by the 
coefficients of the explicit In/x^ terms in (^. For the other masses in (|^), the difference 
between "star" and "bar" is higher order (combined with explicitly two-loop effects), as is 
the implicit yU dependence of the "bar" coupling constants. After these observations, and 
reinstating the (725 S'l terms, differentiating (|^), we obtain 

/^i(p) = it [6P' + (2iV^ - 3 + 12C^|)p - 8iV, 
-dp + \t)P + 3p + + 2 - loop 



where p = r 

Let the right hand side of (H) be denoted Pp. Because of the —8Nc term in (H), there 
is a critical value of p below which (3p becomes negative. And if the starting value of p is 
below this value, as p decreases the derivative becomes more negative, driving p negatuve 
for some value of p, unless some higher order effect intervenes. 

The first higher order effect is the running of the MS coupling constants, which appear 
as coefficients in (^, and the dependence of the lower bound on the cutoff (maximum 
value of ^). One has to integrate the coupled RG equations for five independent "coupling 
constants", gl^ghdhf^ P- Let t = In ^. 
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(22 - lNj)ts 




i9i = ^nNf+\NM 

it = T6^[(3 + 2iV.)y^ - l2C,tst - Ut - Wit] 

Hp) = it [6p2 + (2iVe - 3 + 12C^||)p - 8iV, 
(9 9l I i9l\^,og^ I 3/ g2+gi \2i 



(9) 



The first three equations are integrated trivially. If we neglect the g2,gi contributions 
to the y running, that equation can also be integrated analytically. But if one runs 
up to high scales, the electroweak gauge couplings become of same order as the QCD 
coupling constant; so we use NDSolve from Mathematica to provide an interpolating 
function solution for y'^ which is substituted into the p equation, which is again integrated 
numerically by NDSolve. 

Before giving any results, we discuss the question of their sensitvity to the choice of 
dimensional regularization scale p. In this approach, the scale sensitvity comes from the 
choice of the Electroweak scale fiQ. This enters in three ways: (i) the starting values of 
the gf,t at if: = (/i = po). (ii) the connection between tmax and the nominal value of 
the cutoff, Pmax = A. (hi) the conversion of the found critical Pc(0) back to the ratio of 
pole masses, (^). It is clear that there will be substantial cancellation between these. 

We take rrit = 175 GeV. For the low scale cutoff, we take A = 1 TeV. For a first 
orientation, we keep only the large coupling constants , gs,y, in equations (H), {g2,gi 
0). To check the sensitivity to the choice of starting Electroweak scale po, we solve the 
remaining three equations from @ starting from /io = Mz (tmax = 2.4), and again, 
starting from /io = rrit (tmax = 1-74). We also require as input the initial values glit- 
For muQ = Mz, We take g"^ = 1.483 {as{Mz) = .118). For the Yukawa coupling constant 
we have 



v*{Mw) = {V2Gf)-^ (1 - Ar*) 



251GeV 
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y*^(M^)=2-^ = .972 

V 

Then we use equation (59) of to convert from star to bar, giving the initial value 
y^{Mz) = .970. With these initial values, the critical initial value of p(0) for which p{t) 
falls through zero at t = 2.4 is Pc(0) = .26685 = {M/rh)^. Converting back to a ratio of 
pole masses by gives the critical value M* = 74 GeV. 

For initial scale /iq = rrit, we need glimto) and y'^ijnt). For the QCD coupling constant 
, to our level of accuracy, we can simply run g1 from /i = Mz to fi = rrit using the first 
equation of (|^). This gives g1{mt) = 1.355 {asimt) = .1078). The determination of 
the Yukawa coupling constant is more complicated. We can proceed in two different 
orders. First, we can just run y'^{Mz) to y'^{mt) using the fourth equation of (||) and 
obtain y^{mt) = .905. Alternatively, we can run v*^{Mz) to v*^{mt) (see equation (40) of 
(HI)) convert to y*^ {nit) and then to y'^{nit) which gives y'^{rat) = .877. So we carry 
through the calculation for three different initial values: y'^{mt) = .905, .891(af (7), .877. 
The resulting Mq are (resp) 74.7,74.1,73.5. These results give us some confidence that 
the calculation is not sensitive to the choice of initial weak scale in the range from Mz to 
rrit, as long as the scale dependence of the input parameters is handled consistently. 

Having checked this point, we reinstate the gauge coupling constants in(^. We in- 
tegrate this set of equations, starting at fiQ = rrit , with the additional input gli'fnt) = 
.4239, gf{mt) = .1260. Then the critical initial value of p(0) for which p{t) falls through 
zero at t = 1.74 is .182. Converting back to a ratio of pole masses by equation (|^) gives 
M* = 72 GeV. We take this value as our best estimate of the smallest Higgs mass for 
which the MS perturbation theory is nontachyonic up to scale p equal to one TeV (for 

= 175 GeV). 

Taking the large scale to be the Planck scale {fimax — 10^^, tmaa; — 39) and using the 
same input at fiQ = rrit, we find the critical value of the Higgs mass to be 140 GeV. 

In conclusion, we find that the MS perturbation theory develops tachyonic singular- 
ities (negative mass squared in the MS renormalized propagator) when the dimensional 
regularization scale factor /i exceeds some assigned cutoff value unless the physical Higgs 
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mass exceeds some minimum value, which depends on the cutoff value. 

The requirement that MS perturbation theory not develop tachyonic singularites be- 
low some prescribed cutoff scale is in principle different from and independent of the re- 
quirement that the effective potential not develop a stable minimum at some value of the 
vev much greater than the weak scale. However, until the recent papers of Casas,Espinosa, 
and Quiros 0] it has been generally taken that the requirement on the effective poten- 
tial was practically equiivalent to the requirement that the MS running quartic coupling 
constant stay positve below the cutoff scale. And in the minimal Standard Electroweak 
Model, in conventional renormalization schemes, the condition of positive ratio of squared 
masses implies, and is implied by, the condition of positive quartic and squared Yukawa 
coupling constants (^. In this case the numerical results are not different (up to dif- 
ferences in handling input parameters,etc); and our numerical results are the same as 
those of Altarelli and Isadori But CEQ have argued that consideration of the scale 
dependence involved in the minimization of the perturbative effective potential leads to 
the requirement thatAg// stay positive, where Ag// is not the same as A. Then the numer- 
ical equivalence is broken and the results of the two approaches are different. It is also 
possible that the singular behavior we have found is just a pathology of MS perturbation 
theory, which is interesting in itself, but not justifcation for a restriction on the physical 
Higgs mass. 

It is clearly desirable to have a large scale lattice simulation study of the combined 
Higgs-heavy quark-QCD sector. (Contributions from light quarks and electroweak gauge 
bosons are small, particularly if one doesn't run up to some very high scale). We note 
that a quenched approximation simulation is not adequate for this problem. The term 
in (IP which triggers the possible instability is the —8Nc, clearly a contribution from an 
internal fermion closed loop. 
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